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Abstract 

The main objective of this paper is to obtain an operator realization for the bosonization 
of fermions in 1 + 1 dimensions, at finite, non-zero temperature T. This is achieved in the 
framework of the real time formalism of Thermofield Dynamics. Formally the results parallel 
those of the T = case. The well known two-dimensional Fermion-Boson correspondences 
at zero temperature are shown to hold also at finite temperature. In order to emphasize the 
usefulness of the operator realization for handling a large class of two-dimensional quantum 
field-theoretic problems, we contrast this global approach with the cumbersome calculation 
of the fermion-current two-point function in the imaginary-time formalism and real time 
formalisms. The calculations also illustrate the very different ways in which the transmutation 
from Fermi-Dirac to Bose-Einstein statistics is realized. 



1 



1 Introduction 



The bosonization of fermions has proven in the past to be a very useful technique for solving 
quantum field theoretic models in 1+1 dimensions pQ. Recently this Fermion-Boson mapping has 
been discussed for the case of a finite Temperature, using the imaginary time formalism [2], and 
the equivalence of the massive Thirring model and the sine-Gordon theory was shown also to hold 
at finite temperature. The main purpose of the present paper is to fill a gap in the literature, by 
considering the operator formulation of two-dimensional bosonization at finite temperature. To 
this end we shall use a real time formalism of Thermofield Dynamics |Sl lU El E| • This will enable 
us to give a unified treatment of whole class of two-dimensional models at finite temperature. 

The paper is organized as follows: We begin in section 2 by computing the two-point function 
of the current of free fermions within the imaginary time and real time formalisms. In this com- 
putation we show how the transmutation from Fermi-Dirac statistics to Bose-Einstein statistics is 
realized in quite different ways. These calculations also serve to emphasize later the simplifying 
aspects of the operator thermofield bosonization. 

In section 3 we consider the two-dimensional free massless scalar thermofield. The correspond- 
ing two-point function exhibits two infrared singularities, one similar to the zero temperature case 
and a new temperature dependent one, implying an indefinite-metric Hilbert space on which the 
massless scalar thermofield acts. We then show that one can nevertheless define positive definite 
Wick ordered exponentials of a massless scalar thermofield, provided we associate with them a 
conserved charge (superselection rule). 

In section 4 we consider the two-dimensional Fermi thermofield and compute the corresponding 
two-point function. In the zero temperature limit we recover the known results. 

In section 5 we then use the Wick-ordered exponentials of the free massless scalar thermofield 
as building blocks for the operator thermofield bosonization of the free massless Fermi field and 
show that the Fermi thermofield satisfies the correct statistics. We also compute the fermion ther- 
mofield current by a point-splitting limit and verify that the current satisfies the usual bosonization 
correspondence. As a matter of fact, the well known Fermion-Boson correspondences at zero tem- 
perature 

^irfd^ - l ~ -.d^dy:, (1.1) 

$-fil) -> --^e""3„0, (1.2) 
Miptp -»• - ^— : cos : , (1.3) 

7T 

are shown to hold also at finite temperature. 

In section 6 we illustrate the use of thermofield bosonization by solving the massless Thirring 
model at finite temperature and conclude with comments and outlook in section 7. We leave to 
the appendix the computation of the two-point function of the free massive scalar thermofield, and 
show that the massless scalar thermofield two-point function is recovered in the zero mass limit. 
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2 The Transmutation of Fermi- Dirac to Bose-Einstein Statis- 



tics 

In this section we illustrate how the transmutation from Fermi-Dirac statistics to Bose-Einstein 
statistics is realized in the imaginary and real time formalisms in quite different ways. These calcu- 
lations also serve to emphasize later the simplifying aspects of the operator thermofield bosoniza- 
tion. 

2.1 Current 2-point function in the imaginary time formalism 

In the following it is convenient to work in light-cone coordinates x± — X\ ± ix 2 . Defining corre- 
spondingly j± = ji ±ij 2 , we have for the euclidean two-point current correlator at zero temperature 



<3+(x)j+(y)> = -tr{"f + iS E {x-y)-f + iS E {y-x)) 

d 2 k 



(2*) 



where the overall minus sign arises from the Fermi statistics, and where, after performing the trace 
using 4 = \(l+q~ + 7-9+) and 7^ = r/i = 0, 

II<2(fc) = A J -JL[(p - fc)i + i(p - k) 2 ][ Pl + ip 2 ]D E {p - k)D E {p) (2.4) 



with 

1 



(27T) 

D E {q) 



q 2 (qi + iq 2 )(qi-iq2) 
At finite temperature (|2.4|) is replaced by 

< j+{x)j + {y) >p= £ i / ^^^^^e^^-^f^,^) 

where ue are the Matsubara frequencies U£ = with (3 = 1/T the inverse of the absolute 
temperature, and 

1 f dpi 



flS(A; 1 ,^)=4V - / ^i- 



i(m+|)f]bi-^i-i(m-^ + |)f] 



The evaluation of this expression requires regularization. Indeed, as is well known, formally the 
result depends on the order in which the p\ integration and summation over m is being done. 

Perform first the summation over m 

It is convenient to rewrite this expression as 

2 



z/3 \ ^ 1 r dpi 1 



/3 

3 



Making use of the familiar formula [B] 



77 77 777 77 7 = (tanh7nr — tanh iry) 

m [i m + 2) + lx \ K m + 1) + % y\ x ~y 



we obtain 



Now 



tSDh^q 1 = e(q 1 )(l-2N F (\q 1 \)) 



where Np^qil) is the Fermi-Dirac distribution 

eflftl + 1 ' 



N F (\qi\)= * , : • (2.5) 



Proceeding with the pi integration, one finds 

A similar result is obtained in this way in the T = case [7j. 

Perform first the integral over pi 

Performing first the integral in pi, one obtains 

dpi 1 

2^ [ Pl - i{m + [ Pl -k 1 -i(m-l + \)lf} = 
= -i [6 (2m + l)0(-2m + 2£ - 1) - 0(-2m - l)0(2m - It + 1)] 

Summing now over m one finds, 



i 



it k\ — i « 

This is again the analogue of the corresponding T = result [Zj. 
Averaging now the results we obtain 
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U. ++ (ki,ut) 
As the final step it remains to compute 



(2.6) 



>e p 



h-i 2 -f 



We may write this as 



where 



< j+{x)j+(y) >p= --d\Df{x - y) 

7T 



Df{z) 



(3 J 2n {k\+ 2 -fy 



The result, valid for \z^\ < (3 is well known. Replacing the sum by an integral by using a type of 
Sommerfeld- Watson transform, one finds, 



d 2 k 



Ak-. 



-- t 27c5(k 2 )N B (\k 1 \) 



where NB(\ki\) is now the Bose-Einstein distribution function 

1 



MM 



(2.7) 



e/9|*il - 1 ' 

Continuing to Minkowski space this is just the finite-temperature scalar two point function of the 
real-time formalism. 



Df{z)^iDf(z)=ij 



d 2 k 



-ik-z 



k 2 + 



16 



+ 2m8{k 2 )N B {\k l \) 



so that 



1 



< Tj + {x)j + {y) > p = -cP + iDf(x - y) 



7T 



(2.8) 



(2.9) 



We have thus verified the correspondence (J1.2)) in the imaginary-time formalism, witnessing 
the metamorphosis from the Fermi-Dirac distribution describing the intermediate fermionic states 
to be Bose-Einstein distribution in the equivalent bosonic description! The demonstration for the 

component of the current correlator proceeds in the same way. Since the H — components 

vanish because 7^. = 0, this proves the full equivalence. 

We now show that in a quite different, and surprising way, this metamorphosis is realized in 
the real-time formalism. 
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2.2 Current 2-point function in the real time formalism 

Using the standard rules of the real-time formalism, we have for the ++ component of the current 
two-point function, 

< Tj ++ (x)j ++ (y) >a= -tr (pf + iSf\x - y)^ + iS { p\y - xj) , 

where 

Sf{z) = J d 2 q4Df{q)e- iq - z 

with 

D { P{q) = + 2 m 5(q 2 )N F (\q\) , (2.10) 

q A + le 

where Np(\q^) is the Fermi-Dirac distribution (j2.5j) . Here we use x ± = x°±x 1 and correspondingly 
ry± _ _|_ Noting again that 4 — |(7+<?- + 1-1+) an d 7+ = 7- = 0, we find 

< Tj + (x)j + (y) >p = J ^n ++ (k)e- ik <*^ (2.11) 

with 

n ++ (k) = 4 / ^- 2 {p - k) +P+ Df{p - k)Df{p) . (2.12) 

We next compute explicitly the integral in 1)2.12)1 . Separating into terms independent of the 
temperature (IT++), linear in the Bose-Einstein distribution (II++), and quadratic in iVi?(|g|) (II++), 
we have 

n ++ (A;) = n®(k) + n&(*) + n&(*) 

we have 



n£l(*) 



<i 2 p (p — k) + p + 

(2tc) 2 [{p - k) 2 + ie}{p 2 + te) 



+ - 2 , * (2?rz)iVF(|pi - fci|)g((p - fc) 2 )] 



p 2 + ze 

can be calculated in a variety of ways. A regularization consistent with the conservation 
law d^j^ = yields 



We now make use of the Plemelj's decomposition 

ol _ / 1 



V 



a(x — x ) + ie \x — xo 
in order to rewrite n+]_(/c) and IT++(A;) in the form 



— ine(a)5(x — x ) 



+V (J-^ (2m)N F (\ Pl - hDip - k) + 5((p - k) 2 )} 



d 2 p 



(2m) [e(p + )<f(p_)e(p + - k + )6((p - k)-)[N F {\ Pl \) + N F (\ Pl - h\ 



and 



j2 

n { ? + (k) = 4 / j^(2ni)N F (\ Pl \)e(p + )S(p_)(2ni)N F (\ Pl - k 1 \)e{p + - k + )5{{p-k)^) 



Now, 



V 



j-^y2m)N F (\ Pl \)p + 5(p 2 ) = -±N F (^p) e(p + )6(p-) 



1 



1 



V\ — \(27ri)N F (\p 1 -k 1 \)(p-k) + 5((p-k) 2 ) = —N F 



'\ip-k). 



e(p + -k + )5((p-k)-). 



Since the left hand side of the above expressions are odd functions of p + and p + — k + , respectively, 
they do not contribute to the integral above. In terms of light-cone coordinates we have d 2 p = 
\dp + dp-. Performing the p_ integration we are left with 



n ++ (*) 



i^- W M/^(P + -Me(P + ) 



We now make the change of variable ^ = q + | with k = k + /2. Let us suppose for the moment 
that k + is positive. We may then split the integral in the following way: 



n ++ (A;) = - l - kJ - 



IT k 2 

+ 25{k_) 



— CO 

oo 



dq 



+ 25(k„) I dq 



+ 25{k^ 



dq 



N F {~q -\) + N F {-q + |) - 2N F {-q - \)N F {-q + 
N F (q + k -)+ N F (q - |) - 2N F (q + ^)N F (q - ±) 
N F (q + |) + iV F (| - - 2iV F (g + \)N F (± - q) 



Making now explicit use of the Fermi-Dirac distribution (|2.5|) . we obtain from here 



U ++ (k) 



tt k 2 



+ A5{k_ 



dq 



cosh (3 



dq- 



cosh (3q 



cosh /3| + cosh (3q Jo cosh/5^ + cosh 



which can be rewritten in the form 



n ++ (A;) 



i k+k A 



7T 



k 2 



+ 2S(k. 



I{k) cosh (3— — k 



(2.14) 



where 



/OO 
dq 
-oo 



k 



oo cosh + cosh / 



sinh/3§ 

Putting things together we thus have 



Our analysis was done for A; + > 0. Since our original expression ()2.14|) is symmetric under k — > —k, 
we find that in general 



n ++W = --^UT + 2|fc + |<y(fc_)- 1E - J - 



Noting that 

|Jfe + |<y(ife_) = k + k + 5(k 2 ) 

we finally conclude that 

n++(A;) = ~k+k+ (± + 2Tn5{k 2 )N B {\k l \)^ 

where iYe(|A;i|) is now the Bose-Einstein distribution (|2.7p . A corresponding result is found for 
II (k), while again n + _ = II_ + = 0. We thus conclude that 

< Tj^x) Ju (y) >p= -~d,~d u iD { g\x - y) 

71 

in accordance with the correspondence (jl.2j) for finite temperature. Notice that unlike in the case 
of the imaginary-time formalism, we have now witnessed the actual conversion of the Fermi-Dirac 
to the Bose-Einstein distribution through an integration process. 
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3 Free Massless Scalar Thermofield 



As the first step we shall consider the massless scalar field since it is the building block in the 
bosonization of two-dimensional quantum field theory models. The construction of a quantum 
field theory at finite temperature requires doubling the number of fields degrees of freedom [31 
IH IH1 • This is achieved by introducing fictitious "tilde" operators corresponding to each of the 
operators describing the system considered. This fictitious system is an identical copy of the 
original system under consideration, except for the opposite norm of the corresponding field. This 
entails a doubling of the Hilbert space. 

To begin with, let us introduce the free scalar thermal doublet 1 , 



corresponding to the total Lagrangian density, 

C T = C-C = l - d^d^ - ~ d^4> , (3.2) 
with the corresponding equations of motion (<9± = d ± di), 

d + <9_ <f>(x) = , <9+<9_0(x) = . (3.3) 

In view of the equations of motion ()3.3)1 . the two-dimensional free massless scalar field can be 
decomposed in terms of left- and right-movers (x = x° ± a; 1 ), 

<P(x) = <P L (x + ) + <P R (x-), (3.4) 

and similarly for (j){x). In order to simplify the notation, we shall omit the subscripts, which are 
taken to be implied by the arguments x . 

In the quantized theory, the two-dimensional Boson field at zero temperature is described by 
the field operator , 

(dp) (f p {x) a{p l ) + f^a^ip 1 )). (3.5) 

-co N ' 

f p (x) = , (3.6) 

(dp) = ^ . (3.7) 
/ (2tt)(2|^|) 



The tilde conjugation is defined by the property 

(ca) = c*a , (3.8) 

1 The conventions used are: 

7°=(J 1)^=^ 1 Q ),l 5 =l 7 1 ,e^ = l,g 00 = l,x ± =x°±x\d ± = d ±d 1 . 
For the free massless scalar field <f>(x) = <f>(x~) + (j>(x + ), and for the pseudo-scalar field <p{x) — <p(x~) — (p(x + ). 
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such that the tilde conjugated field is given by, 

(dp)(r p (x)a( P 1 ) + /,(x)aV))- 



(x) , <f>{y)} = 



At T = these fields are independent, 



The Fock vacuum state is, 



|0,0) = |0) <g> |0) 



with the property, 



a(p 1 )|0) = , a(p x ) |0) = 0. 
Decomposing the free massless fields into left- and right-components, we obtain 



dp 1 




f ai-p 1 




\ aip 1 ) 


dp 1 


>;(^) 


fdi-p 1 


J Ait Ip 1 ! 


{ dip 1 ) 



d^p 1 ) 



where, 



and 



[0(x ± ),0(^)] = [0(a; ± ),0(^)] = O. 

Note that in Eqs. (jH.lHj) and 1)3.14)1 the limits of integrations are [0, oo). 

In Thermofield Dynamics the temperature-dependent vacuum is defined by, 

|O(/?)) = E/b(0b)|O,6), 
where the unitary operator Ub(0b) is given by, 

U B (6 B ) = e~ iQ{eB) = e~ /- + r dpl (^ pl ) a ( pl )- at ( pl ^ t ( pl )) 9s(|pl| ' /3) 
and the Bogoliubov parameter 9b{\p 1 \,I3) is implicitly defined by 

sinh^dp 1 !;/?) 
cosh 6 B (\p x \;P) 



1 



with the Bose-Einstein statistical weight is given by, 
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N B (\p 1 \;P)=smh 2 6 B (\p 1 \;P) = ^ _ ^ . (3.21) 
The correspondingly transformed annihilation operators are given by, 

atf;0) = U B (-e B )a(p 1 )U B (6 B ) = a^ 1 ) cosh^dp 1 ]; /?) - ^(p 1 ) sinh^dp 1 ]; /3) , (3.22) 

atf;P) = U B (-6 B )a{p l )U B (6 B ) = a{p l ) cosh 6 B {\p l \; (3) - a"* (p 1 ) sinh B 9(^1; (3) . (3.23) 
The vacuum state at finite temperature satisfies 

aip 1 ; (3) |0(/3)> = , a{p l - (3) |0(/3)> = . (3.24) 

Here, and in what follows, p is defined to be p = {p 1 \ . 

The thermofield operators, that act on the Fock vacuum |0,0) are given by, 

<P{x ± -(3) = <p {+ \x ± -(3) + <p { ~\x ± -(3), (3.25) 

${x±;P) = ft + \x ± -(3) + ft-\x ± ] (3), (3.26) 

with 

^\x±-(3) = - 4>f\x ± ] (3) , (3.27) 

^\x±;(3) = 4>f\x±-(3) - ct>f\x^(3) , (3.28) 

<Pi ± \x ± ;(3) = f(dp) (1^)^)) ™ she B(p;(3), (3.29) 

= f(dp) (f^a^p)) sinh ^^^)' ( 3 - 3 °) 

Jf=V;/3) = jT(*) sinh^(p ;/3) , (3.31) 

^V;fl = jf (*) cosh M ^). (3.32) 
The commutators of the thermofield are the same as that at zero temperature, 

[00^; (3), 0(y ± ; = [0 c (x ± ; 0^; + [0 i (a; ± ; /3), s (y ± ; = [00^), 0(y±)] , (3.33) 
so that 

[0(^5 (3), 0(^5 = [0(a: ± ), 0(/)] , (3.34) 
11 



and 



[^x ± ;P),<P(y ± ;P)} = 0, (3.35) 

As in the zero temperature case, the massless scalar theme-field carries an identically conserved 
topological current associated with the Gauss' law 

J' i (x]/3)=d v €^ ( p(x]/3), (3.36) 
and the field (p(x; (3) is the dual of <j)(x; (3), 

e^d uV {x;(3) = d^<f){x;(3). (3.37) 
The corresponding charge is formally given by 

Qp = J J°{x;f3)dx 1 . (3.38) 
3.1 Scalar Thermofield Two-point Function 

For free fields the n-point functions are determined by the two-point function. The diagonal 
contribution to the Schwinger two-point function is given by 

<0(/?) 10(^)0^) |0(/?)> = 

where in our simplified notation, 
Di + \x ± - y^, P) 

Di + \x ± -y ± ;(3) 

We get, 

(O,O|0(a; ± ;/5)0(y ± ;/3)|O,O)= J DW( a ; ± -y ± ) + ^ f° ^ [cosp^ - y±)] N B (p; (3), (3.44) 

Z7T JO p 

where the zero temperature two-point function is given by, 

D M( X ± - y ±) = - _L enMx* -y ± - ie)] . (3.45) 



(O,O|0(rn ± ;/3)0(2/ ± ;/5)|O,O> (3.39) 
[^\x ± -p) ) ^-\y ± U3)} + [^(x ± ;f3)Ji-\y ± ;f3)} (3.40) 
D i+)(a± -y±;(3) + Di + \x ± - y±; (3) , (3.41) 



r fp(x±) ^ c ° sh2 ^ ® ' (3 - 42) 
r WW f ? {x±) fp(y±) sinh2 Mp; p) ■ (3 - 43) 
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where // is an arbitrary infrared (IR) regulator. Besides the usual infrared divergence, the integral 
in ()3.44|) introduces an additional divergence. The integral appearing in ()3.44|) is given by, 2 



^- [ — [cosp(x ± -y*)] N B {p;P) 
In Jo p 



-M^-V*) - -Inlsinh^-^j - UJ- + — f° . (3.46) 



4tt ' An \ (3 j An n In Jo p (eP p — 1) 

The integral appearing in (J3.46|) is infrared divergent. We shall introduce an infrared cut-off \j! 
and define 

which corresponds to the mean number of particles having momenta in the range \jjf, oo], 

z(/3,//) = (0;P\JV B \0(P)), (3.48) 



-Vb = / -a\p)a(p), (3.49) 



where 



V P 

is the operator for the number of particles. For a fixed temperature the asymptotic behavior of 
the integral (|3.47j) . besides the logarithmic singularity, exhibits also an algebraic singularity, 

z(fi'^0;(3)^^ + Un((3ii'). (3.50) 

As we shall see, the integral ()3.47|) plays a crucial role in determining the selection rule for the 
correlation functions of Wick exponentials of the free massless scalar thermofield. 

From (|3.44|) . (|3.46|) and (|3.47|) the two-point function for the left- and right-movers seen to be 
given by , 

(O^I^/M^/mO) = -^ln{^sinh ^ (x± ~^~* £) } + ^*(&a0- (3-51) 
For the field <j)(x ± ), we obtain 



2 Gradshteyn - 4 Edition, pg. 495, 19 m = 1; 

00 cos(aa;) - cos(fcc) gLt 1 a sinh ^ 



o e#*-l x 2 6sinh^ 
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The limit e — > + was introduced as a prescription to properly define the zero-temperature limit 
of the above expression as a tempered distribution. For the space-time contribution, the short- 
distance limit x ~ y coincides with the T — > limit, 

(0,010(0^; ^0(^5 ^16,0) - -i ;: ln{i/i(ar t -3/ ± -ie)}. (3.53) 
The two-point function for the scalar thermofield 4>(x; (3) = <ft(x + ] (3) + <p(x~; 0) is given by, 

(O,O|0(x;/3)0(y;/3)|O,O) = (3.54) 

7T 



— In < — u 2 f — ) sinh 

47T 



/3\2 . r7r(x + — y + — z'e) 



sinh 



/3 /3 



7r(x — y —it) 



In the zero-temperature limit we obtain the two-point function of the free massless scalar field, 

(O,O|0(x;/?)0(y;/?)|O,O> T ^o - (O|0(x)0(y)|O) o = - ^ In { - f (x-y) 2 + *e(x -y )} . (3.55) 

The fact that the short- distance behavior of the thermal two-point function is the same as that 
for zero temperature is a crucial factor for the success of the thermofield bosonization scheme. 
The appearance of the additive temperature dependent arbitrary constant in (|3.51|) . z, is of no 
immediate significance since the correlators of the massless scalar field in two dimensions does not 
satisfy positivity. As we shall show in the following section the dependency on this constant of the 
Wick ordered exponentials will disappear upon imposing suitable superselection rules. 

The off-diagonal contribution for the two-point function is given by 3 , 

±. a\Zt..±. flMn n\ 1 f°° _ ^± ..±1 M - m ^ d P 



(O,O|0(x ± ;/3)0( 2 / ± ;/3)|O,O)= -— / cosp(x ± - y±) e~ N B (0,p) ^ 

2n Jo p 

1 f°° cosp(x ± — y^) dp 



47r jo sinh |p p 

= l l * cosh n{x± p y±) ~h mK) > (3 - 56) 

where f(/3, k) is the integral, 

m , K)= r i *= r± NMp;0)e ¥. (3.57) 

Jk smhfp p Jk p 

where k is a new infrared cutoff. This integral also play an important role in the selection rule of 
the off-diagonal correlation functions of thermal Wick exponentials. In the limit T — > or x ~ y, 
we get for the space-time dependent contribution, 

(0, 0100^; PMy*; (3) |0, 0) T ^ - , (3.58) 

as expected. 

3 Gradshtcin pag.494, eq. 18. 
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3.2 Wick Ordered Exponential and Selection Rules 

Let us consider the Wick ordered exponential of the free massless fields cjy^x^) and (f)(x ± ), 

W(x; A) =: e iX ^ x) := e iA0( ~ )(x) e iA0(+)(x) , (3.59) 

W{x; A) =: e~ iX ^ x) := e -^*(-)(*) e -iA#+i(«) _ ^ 6Q ) 

The Wick exponential has to be understood as a formal series of Wick-ordered powers of the field 
at the exponent such that, 

(0\W(x; A)|0) = (0\W(x; A)|0) = 1 . (3.61) 
The unitary operator U B (6 B ) ac ts on the creation and annihilation components as 

U B {B B ) <t> {± \x) U B \6 B ) = 4>f\x; 0) - ^\x; P) , (3.62) 

U B (6 B ) fi+\x) U B \9 B ) = 4>f\x- P) - ^\x; (3) . (3.63) 
The thermofield Wick exponentials are then given by, 

W(x; (3, \) = U B (-9 B )W(x; \)U B (6 B ) = (V A< ^~W) a^+W) ) (^ix^^-ix^^p) ) 

= Z(f3,n',\ 2 ) ( e *^£~W)-**$~W) ^Lix^&fi-ix^ixtf)^ ^ 

W(x;P,\) = U b (-6 b )W(x;\)Ub(6b) = (V* A ^~W)+^<+W)) ^-i^^+ix^x-^ 
where 

x2 



Z(/3, / u',A 2 ) = e- A2 ^ +) (^'^" ) ^|^ o = e -^f~% NB{k ' P) , (3.64) 



that is 



Z(/3,/i',A 2 ) = e~£ z{M , (3.65) 

where z((3,fi') is the integral appearing in the two point function (|3.46|) . The thermal Wick 
exponential in then given by, 

W{x- (3, A) = Z{P, //, A 2 ) : e iA ^ ;/3) : , (3.66) 

W(x;P,X) = Z{P,fi',X 2 ) :e - iA ^) : . (3.67) 
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Thus, we obtain, 



W)\W(x;\)\0(P)) = (0,0|W(x;/3,A)|0,0) = Z((3, fi! , X') 



(3.68) 



(0{(3)\W(x;X)\0(f3)) = (0,0 \W(x; (3, A) | 0, 0) = Z((3,fi', A 2 ) . (3.69) 

As we shall see, the factor Z(f3, fi', A 2 ) plays a role of a wave function renormalization for the 
bosonized Fermi thermofield. The exponential is IR dependent on the parameters fi and fi', such 
that 



■ ' e i\(f>(x;py 













(3.70) 



where iV^' {N v y) stands for "normal ordering" with respect to the masses fi, fi! [y,v'\ 

In order to display the role played by the cut-off dependent term Z([3,fi', A 2 ), let us consider 
the correlation functions at finite temperature of the Wick exponential, 



Wix^Xj) = : e 4A ^ } : 
The diagonal contribution is given by, 



(3.71) 



J] W{ Xj ;\ s )\O(0)) = (0,0|n tn*;;/3,A;)|0,0> 

3=1 3=1 



e -£*{M ZU a 2 



-/ ■ ■ "/ J] e 

i<j 



e -±z{[3,fi') E"=iA 2 -Q g-AiAHDi+'^-^^+DW^-^;/?)} 



Using the expression (J3.51JI for the thermal two-point function, we obtain, 



(3.72) 



w)! n w ( x 3^j) w)) = e ^0(^=1 Aj) x 

i=i 

t sinh 7r (^~^-^) ] ^ (3 73) 

In order to ensure that the correlation functions of the Wick exponential at finite temperature be 
independent of fi', we must impose the selection rule, 

n 

]TAi = 0. (3.74) 

i=l 

Thus, we can write, 
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w)in w{ Xj] \)\m) = <w (0,01 n wixthXj)^) 

3=1 3=1 



Itt 



7T /? 



Xj - n -/? vr(xi -x i -ie) A,A ' 



- <w„ n (a*) 8w n [£ j ' } 47r ( 3 - 75 ) 

7=1 i<3 P 

As in the T = case we may thus define a multiplicatively, infrared renormalized Wick exponential, 

Wife A, ) = (fi^WfaXj), (3.76) 

such that, 

WnixfifrXj) = (ii)*Z(p,S, A?) : e iA '«**0 : . (3.77) 

It is interesting that one and the same selection rule J2iXi = allows us to eliminate the dependence 
on the infrared regulators fi and fi' . Thus, provided we associate with the Wick exponential a 
conserved charge A, the correlation functions of the thermal Wick exponentials (J3.77J) are mapped 
into Wightman functions belonging to a positive metric Hilbert space. In a similar way for the 
correlation functions of the Wick exponential W{xj \ [3, Xj), 

A 2 

Wnix^frXj) = (fi)^ Z{(3,fj!,X)) : e'^^P) . . ( 3 . 78 ) 

Using Eq. (|3.51j) . the correlation functions of the Wick exponential W(x; f3, X) can be obtained 
from (|3.75j) by making ie — * — ie and introducing a phase factor 



En ^i^j 
*<i ~ . (3.79) 

The off-diagonal correlation functions are given by, 



<o(0)i n wto^) fi wiy^x^om =^{& x >) +te m \) } 

i=l i=l 

e & { (e: a,) 2 + (e; a,) 2 } e - ± m k) (e? a, ) ( a 



X 



TT [ — sinh 

;J7 L 7T 



i<fc 







4-7T 



TT I" — sinh 



Aj A^ 

4-7T 



it,. no 

IT [ cosn 



'•7 



/ \ AjA ? - 



(3.80) 
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Again, the selection rules that makes the off-diagonal correlation function simultaneously indepen- 
dent of /i, // and k are, 



E ^ = » E = o 



(3.81) 



i=l 



and we get for the IR renormalized Wick exponentials, 

n m 



i=l 



rr f • ^ • u ^(^ _ 



J][ \i — sinh 



/3 



AjAfc JTi 
47T 



^ <Vi-Vt) 



A, A 



TT [ z — smh 



77 /3 



4-7T 



X 



It, fib 

] { [ cosh 



I -J 







In 



(3.82) 



4 Free Massless Fermi Thermofield 

Let us introduce the fermion doublet jlj 

/ ^ 



(4.1) 



and the corresponding total Lagrangian density, 

C T = iV-fdfl = £-£ = ilhM - ( -i^-fd^ . (4.2) 
The free massless Fermi fields ip and ip are given in terms of right and left spinor components, 



ip(x) 



ip(x) 



ip(x~ 

ip(x^ 
ip(x~ 



where 



JO \/27T P 



i/j(x^ 



V27T 

dp 



The unitary operator taking one to the thermofields is now given by, 



U F (8 F ) = e 



J™ co d P e F (\p 1 \,f3)(b(p 1 )b(p 1 )-bH P 1 )bHp 1 ) + dip 1 )d(p)-dH P 1 ) £ (p 1 ) 



(4.3) 
(4.4) 

(4.5) 
(4.6) 

(4.7) 



and the corresponding transformed annihilation operators are 

fcfa 1 ;/?) = b{p l ) coa9 F (p',P) - b\p l ) mn9 F (p;P) , (4.8) 

b{p\{3) = b{p l ) cos F 0(p; 0) + 6+Cp 1 ) sin0 F (p;/3) , (4.9) 

with similar expressions for dip 1 ; (3) and dip 1 ; (3), where the Bogoliubov parameter 9 F ip,P) is now 
implicitly defined by, 

cos9 F ip;P) = 1 (4.10) 
V 1 + e~P p 

e -/3p/2 

BnOrtoP) = , (4.11) 
V 1 + e~ pp 

and the Fermi-Dirac statistical weight is given by, 

iV F (p;/3) = sin 2 ^(p;/?) = — . (4.12) 
The fermion thermofields are then 



^(^ ± ;/5) = /°° ^ (&(=R>) cos^(p;/3)-fet( T p)sin^(p; 

JO V27T 1 V 



dp 
/2m 

+ (>(TP) cob6 f (p;P) - %p) sin0 F (p;/?)) }. (4.13) 



and 



*l>{x t ',P) = /°° ^{£(x±)(&( Tp ) cos^(p;/3) + 6t( T p) sin ^(p ; 

JO V27T L V 

+ / p (x ± )(dt( Tp ) cos F (p ;/ 9)+tZ(=Fp) sin ^(p; /?))}. (4.14) 
Besides the canonical anticommutation relations we have, 

Wx; P) , 4>iy; P)} = , {^(x; /?) , /?)} = . (4.15) 
4.1 Fermion Thermofield Two-point Function 

The diagonal contribution for the two-point function of the thermofield fermion is given by, 

(0,0mx ± ;P)^iy ± ;P)\0 } 0}= — ^ - — / dpN F ip; P) sinp^ - y±) . (4.16) 

2 % 7r(x ± — y ± ) lit Jo 

The first term in ()4.16|) is the two-point function at zero temperature, 

(01^)^)10)0 = 7- -4 IT • ( 4 - 17 ) 

2i^ix ± — y ) 
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The integral appearing in the two-point function (j4.16|) is evaluated to be 4 

1 1 



I POO 

— / dp Np(p; (3) sin p(x — y) - 
in Jo 

We obtain after appropriate ie prescription, 

(o,o|^;WV;0)|o,6) 



2iTr(x — y) 2if3 sinh[|(x — y)] 



2i (5 sinh[f {x± -y ± -ie)}' 



(4.18) 



(4.19) 



(4.20) 



2i (3 sinh[|(x ± -y ± + ie)] ' 
The diagonal thermofield fermion correlation function satisfies the asymptotic factorization 



property 



lim (0,0|^(z + A;/?)^;/3)|0,0> - (0, 0|V(x; /3)|0, 0) 

A— >oo 



For T — > as well as x ~ y at fixed T, we obtain, 

(o^^/^V^^o) - 

For the off-diagonal contribution we find 5 , 



1 



2i ixyx^ 1 — y ± — ie 



(4.21) 



(4.22) 



(COlv^^M^^o) 



1 Z" 00 (3 

-/ cosp{x ± -y ± )N F {f3,p)e-^ p dp (4.23) 

7T JO 

1 r°° cosp(x ± — y^) 
2ix Jo cos |p 



dp 



2/3cosh[f(a; ± - y±)] ' 



In the limit T — > we get, 



(0,0|^(x ± ;/3)^( 2 / ± ;^)|0,0) r _ - 



(4.24) 



as expected. The off-diagonal two-point function also satisfy the cluster decomposition property. 
Moreover, 



(0,0mx ± ;(3)^(y ± ;f3)\0,0) = 0. 



(4.25) 



4 Gradshtein - 4 Edition, pag. 481- eq. 3.9911-1, 

sin ax 1 



dx ■ 



>J 3x + l 2a 2/3sinh 



an 





5 Gradshtein, pag. 503, last integral. 
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5 Thermofield Bosonization of the Free Massless Fermion 



The scale dimension of the thermofield Wick exponential is d = j-. In analogy to the T = case, 
we seek a representation of the thermal Fermi field in terms of the IR renormalized thermofield 
Wick exponential of canonical scale dimension A = 2i/7r: 

^5/3) = Wr^-P) = CZ c (p,ti') : e 2 ^^) : , (5.1) 
^{x ± ;(3) = W R (x ± ; (3) = CZ c (f3,/x') : e ~ 2l ^ x± ^ : . (5.2) 

where 

W R (x ± ;f3) = W R (x ± ;f3,2^) (5.3) 
Z c (pJ) = Z{fi,n\±K) (5.4) 

and 

Taking into account the super selection rule associated with the Wick exponential of ther- 
mofields, one readily verifies that the set of all two-point functions of the Fermi thermofield, 
obtained in the previous section, are recovered from the bosonized expression (J5.1j) . In particular 
property ()4.25j) follows from the superselection rule for the Wick exponential of the Bose ther- 
mofield. The minus sign appearing in the two-point function for the thermofield ip(x; (3) arises 
from the phase factor ()3.79jl . 

5.1 Statistics 

Let us consider the anticommutation relation between the same spinor components labeled by ±. 
From ()5.1|) we obtain, 

^(x ± ;/3)^(y ± ;/3) = C 2 Z 2 c {P,n') : e 2i ^ x± ^ :: e 2l ^^ y± ^ : (5.6) 

= C 2 Z 2 (/3,^) : e 2i ^[<t>(* ± \P)+4>(y ± -M . e -4*{Di + \x±-v±-,P)+Di + \x±-y±-,l3)} 

= C 2 Z*{{3,ii') :ilj{x ± ;f3)^{y ± ;(3) : Up t ^ ( x ± - y± - ie)\ ) . 
From the antisymmetry of the sinhz under z — > —z it follows that, 

{^(x ± ;/?),^( 2 / ± ;/3)} = 0. (5.7) 

In a similar way, we obtain, 

2i(3 ( sinh[JOr± - y ± - i ej] ~ smli^x* - y ± + it)} ) ' 

The non zero contributions in Eq. ()5.8|) arises from the neighborhood of the point x = y. This 
time we get, 
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{V^; (3) , ^{y ± - (3)} = - y±) . (5.9) 
For the tilde-fields we analogously find, 



{^(x ± ;/3),V(l/ ± ;/3)} = 0, (5.10) 

(3) , ^(yi; £)} = ^(aji - y±) . ( 5 .11) 

In order to ensure the correct anticommutation relations between ip(x + ;/3) and ip(x~;{3), be- 
tween ip(x + ;j3) and ip(x~] (3), as well as the anticummutativity of tilde operators with non-tilde 
ones, Klein factors [HI [TT] must be introduced. To begin with, let us introduce the charges, 



and the Klein factors 



/-boo 
<9 0(^/9) dx l 
-00 

Q± = / d <j>(x ± ;f3)dx 1 



Kg = e 2 e 



Since 



Qt = U B (-6)Q ± U B (6), 
the Klein factor Kg is the transformed Klein factor of the Fermi field at zero temperature, 

K% = u B (-e)K ± u B (e). 

Starting from eqs. (|5.1|) and (|5.2|) the thermofields are redefined as, 

Hx + -f3) = W R (x + ;P), 

^(x-;(3) = K+W R (x-;(3), 
4>{x + -(3) = K+KgW R (x + ;f3), 
${x-;P) = K+Kg-K+W R (x-;(3). 



5.12) 
5.13) 

5.14) 
5.15) 

5.16) 

5.17) 

5.18) 
5.19) 
5.20) 
5.21) 



These fields obey the correct anticommutation relations. The factor Kt in (|5.19|) and Kt in (|5.21j) 
ensures that for different (left , right ) spinor components we obtain, 

{^(x ± ;/3), V(l/ T ;/5)} = 0, (5.22) 
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{^■,(3),^; (3)} = 0. (5.23) 

The factors K% in (f5~20jl and (f5~2T^l ensures the normal anticommutativity between tilde ther- 
mofermion components and non-tilde ones, 

{Hx ± ;P),^P(y ± ;P)} = 0, (5.24) 
{^(x ± ] p),ij(y^-p)} = 0. (5.25) 

5.2 Fermionic Current and e-expansion 

The vector current corresponding to the fermion doublet (j4.1j) is given by, 

Tj7y = J> i -J*, (5.26) 

where, 

jv = , (5.27) 

J^=JyY. (5.28) 
Let us define the finite-temperature left and right current components 

J±{x-(3) = ij\x ± ;P)ij(x ± -p), (5.29) 
JpfaP) = ^(x ± ;(3)^j(x ± ;(3) } (5.30) 

by point-splitting as, 



J ± (x-{3)= lim {^{x ± +e ± ;(3)^j{x ± ](3) - (0,0\^{x ± + e ± ;(3)^j{x ± ;/3)\0,0)} , (5.31) 
and similarly for J^(x] (3). We obtain, 

J ± (x;i3) = ^Z 2 J(3, f j,') lim f : e-^^+^P) . . e 2iV*<K*±ifi . -(y.E.V.)A = 

— 2%(fl,H') lim { : e - 2 *v / ^[^ ± +e ± ;/3)-^ ± ;/3)] . e 4 7rD( +) ( £ ±;/?) + 4 7r5i +) ( e ±;/3) _ (y.E.V.)^ } , 
27T 6 £±^0 ^ ' 

(5.32) 

J±(x;/3) = JL Zl (j3,fJL') lim ( : e ^^(- ± +^) . . e - 2^(^/3) . _ | = 

27T £±^0 ^ > 
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-^Z 2 c (P,fi') lim I :e 2^[^ ± + £ ± ;/3)-^ ± ;/3)] . e 4,#( £ ±;fl + 4, D W( £ ±;ffl _ (m) ) (5 33) 

Since the short-distance behavior of the thermal two-point functions (|3.51j) and (J3.52j) is the same 
as that for zero temperature, and the superselection rule is identically satisfied, we obtain the /i 
and /j! independent result, 

J(x ± ;(3) = lim — 1 —{ :e -2*v^ ± s.±*(* ± ;fl : -(V.E.V.)*}, 

1 ^±0(^/3), (5.34) 



J(x ± ;f3) = lim -^——\ :e ^^d x± 4>{x^p) . _(y.£/y) } 



' d^^x^P), (5.35) 



Defining the pseudoscalar thermofield (p(x; (3), such that, 

e lw dr<p(x;P) = d lt <Kx;P) (5.36) 
the vector currents can be written as, 

J" fa (3) = : tfix; PW{x- p):=-J- e ^ d v <pfa 0) = - ±= & <j>(x; 0) , (5.37) 

'7T V 7T 



^(x; /3) = : ^ (x; /3)7"fe /?) : = - A= ^ d u <p(x; 0) = - ^= d" 0(x; 0) . (5.38) 

V 71 " v 71 " 

From ()5.34|) . the current two-point function is found to be given by, 

(O^ix- 0)^,0)10,0) = - d x± d y ±(O,O\<f>(x;0)(f>(y,0)\O,O) 

71 

(5.39) 



A0 2 wnh 2 [^(x ± -y ± 



In the zero temperature limit we recover the zero temperature correlation function, 

< I^)^)I > = 4^ {3 ±_l±_ i€) 2 ■ ( 5 - 40 ) 

An alternative way to obtain the bosonized form of the current ()5.34j) is to consider the ex- 
pectation value of the T = current with respect to the \0,O) vacuum. To this end, let us again 
define the current two-point function in terms of a point-splitting, 
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<0(/3) I J^x) |0(/3)> = 



lim W)\(^(x ± + e ± )i;(x ± ))(^(y ± + 6 ± )^(y ± ))\0(f3)) = 

(|-) 2 Z*(f3 } yl) £± lim_^ o (0, 0| ( W*(x ± + e ± ;(3) W^; (3) ) ( W*(y ± + 6 ± ;(3) W{y±- (3) ) |0, 0) . 

(5.41) 

Since the superselection rule is identically satisfied, the dependence on fi and // cancels and we 
get, 

(|-) 2 £± lim^ o ( -^L ) (o, 0| ( 1 - 2 % V^e ± d x ± <j>{x^ (3) ) ( 1 - 2 i d y± 0^; (3) ) |6, 0) 

= -(0,0\d x± <P(x ± ;f3)d y± 0(y ± ;f3))\O,O) - * (5.42) 

The singular term appearing in (|5.42j) can be removed by redefining the point-splitting by sub- 
tracting the contribution coming from the vacuum expectation value of the individual current 
operators, i. e., 

(o^ij^j^ioos)) = 

e hm ±Q | (0(/3)| (v>V + e±) ^ (a ±) ) ( ^t (y ± + 5 ±) ^ (y ±) ) | (/3)) - 

- (0C9)| ^ V + £± ) VK**) W)| V> V + 5± ) tfG/*) W)> } • ( 5 - 43 ) 

The same for j^a;*; (3). 

5.3 Bosonized Lagrangian 

The total Lagrangian density of the free massless Dirac thermofield, 

C-C = z^(x;(3h fl d^(x;(3) - ( - i^(x; (3)^d^(x; (3)) , (5.44) 

d^(x; (3)d»(j){x- fyfe (3)d^(x; (3) : . (5.45) 



C T = i*(x;P)>fd li *(x;P) 
is equivalently given by, 

r>bos 

- 2 
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The demonstration follows along the same lines as that given in Appendix E in 1J for the zero 
temperature case. To this end, we must consider the Wilson short distance expansion for the total 
Hamiltonian density 



H T = iV^dtV = H-H = iip^diip - ( - iip^diip) + h.c 



(5.46) 



Since the short-distance behavior of the thermofield is the same as that of zero temperature, we 
find 

H(x; f3) = \ ($(x + e; firfd^x; (3) + h.c.) - {V.E.V.) P 



For 7i we find, 



1 r 



: (d <P(x;(3)) 2 : + : (d^x; (3)f : 



(5.47) 



H(x;j3) = --(^(x + E^fti^d^ixifl + h.c.) - (V.E.V.)/3 



— e 



l 7T 



: (d^{x-(3)f : + : (d^x-13) 



(5.48) 



The phase factor e~ t7T in 1)5.48)1 . that cancels the global minus sign, comes from the phase (JH.79)) 
associated with the Wick exponential W(x;/3). The total Hamiltonian density H,t = TC — Ti. 
corresponds to the total Lagrangian density ()5.45|) for the free massless scalar thermofield. 



6 Massless Thirring Model 

In this section we illustrate the use of the thermofield bosonization by solving the massless Thirring 
model, which is defined by the Lagrangian density 

c T = ifrfdrf + \g (^VX^) - ( - + \g {hy){h^ • (6.1) 

To begin with, let us consider the operator solution of the (zero temperature) massless Thirring 
model, 

^0) = —j= \VA ■ e v ) :, (6.2) 

V 71 " v ' 

fe) = ± (/,) ^ (a2+52) : e~* («^H*«*>) : , ( 6 .3) 

where {(p) is the dual of <fi (4>) and \i is the infrared regulator of the two-dimensional free massless 
scalar field. The /i-dependent factor in front of (|6.2|) is needed in order to obtain in the IR limit 
— > well defined Wick exponentials living in a positive-metric Hilbert space. The corresponding 
thermofield is given by, 
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il)(x; /3) = -= Ifi) Z(P,fi):e\ / : , (6.4) 



where 

Z\fl,n) = e V /, (6.5) 

and similarly for ip(x; 0). We shall consider in particular 

5=-, (6.6) 
a 

in order to obtain in the zero temperature limit a canonical fermion field with Lorentz spin s — |, 

a 5 1 . 

Taking into account that the short- distance behavior of the thermofield two-point function is 
the same as that of zero temperature, we use the Mandelstam prescription ^31 E] and define the 
current through the point-splitting, 



J»(x;P) = Km T( £ )[8l + —8l \ $(x + e; 0) -f^x; 0) - (0, 0| ${x + e; /3) 7 ^(x; /3)|0, 0) 

e=(0,e 1 ) 



where, 



One finds, 



^) = (-/. 2 e 2 ) : 



a 



J»(x;P) = - — e'"'d v <p(x;P) 

Z7T 



J»(x;P) = - — e'"'d v !p(x;P). 

Z7T 



(6i 



(6.9) 



(6.10) 



(6.11) 



The quantum equations of motion are given in terms of the symmetrized form of the interaction 
term of the Thirring current with the Fermi thermofields by, 



i^d^(x;P) = -gfafaPW, ^(x;P)} 

vfd${x; P) = g{Mx; P)^ , P)} , 
with the following identification 

n 47T 



a 



(6.12) 
(6.13) 

(6.14) 



by, 



The unmixed correlation functions for the thermofield spinor components ip a (x; P) are given 
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(0,0|^(x 1 ;/3)---^(x n ;/3)^( 2/i; ^...^t (2/n . /5) |o ) 0)= (27r)- n (z) 2 ^ ^ x 




which represents the generalization of the T = result to T ^ 0. 



7 Conclusion 

Using the framework of Thermofield Dynamics [Hj we have shown how the familiar T = bosoniza- 
tion formulae for fermions in 1+1 dimensions generalize in a natural way to the case of non-zero 
temperature. Despite the appearance of now two infrared regulating parameters it was shown 
that the superselection rule already familiar from the T = case guaranteed the independence 
of physical quantities of both these parameters. The "Mandelstam representation" for thermal 
fermions now allows us to obtain in a compact way the exact and complete solution of a number 
of integrable models involving massless fermions in 1+1 dimensions. 

The next step would now address the problem of massive fermions in 1+1 dimensions. Do 
we again have an equivalence between the Massive Thirring model and the sine-Gordon theory 
? We expect that this is indeed the case, since such an equivalence has been demonstrated for 
correlators of fermion and sine-Gordon thermal fields in the imaginary time formalism [2]. Taking 
a Mandelstam representation in terms of a thermal sine-Gordon field as a starting point, can we 
prove, in a way analogous to what was done in the T = massive Schwinger model J2]> that the 
bosonized thermal Dirac field satisfies the free massive Dirac equation ? These and other questions 
will be subject of future investigations. 
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Appendix 

The purpose of this appendix is to show that the theory of the free massless scalar thermofield 
can be constructed as the zero mass limit of the free massive scalar thermofield theory. To begin 
with, let us consider the free massive scalar thermofield, 

E(x; (3) = U[9 B (/3)\ E(x) U[9 B (J3)] = 
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7-oo VP + m 2 i/p^K ^ cos1i Mp; _ a (P) sinh0 B (p; 

+ /* (x) (a 1 " (p) cosh # B (p; /3) - a (p) sinh 6 B (p; /3) ) } • 
The thermal two-point function is given by, 

(0, 0|E(x; P)E(y; (3) |0, 0) = D^\x -y;m) + l(x - y; /?, m) , 
where — y; m) is the corresponding zero temperature two-point function, 



(A.l) 



(A.2) 



D^ + \x — y;m) — — 



dp 



{< 



(x -y°)y/p 2 +m 2 + (x 1 -y 1 ) 



47T Jo \Jp 2 + m 2 
and the last integral in (|A.2|) can be written as, 



+ e 



j/')^/p 2 +m 2 — (x 1 — y 1 )p 



(A.3) 



1 r 00 
Z(x-y;P,m) = — \ 
In Jo 



dp 



\Jp 2 + m' 2 



N B (w, (3) { cos [(ar° - y°)^Jp 2 + m 2 + (x 1 - y^p 



cos 



(x° — y°)\Jp 2 + m 2 — (x 1 — y 1 )p | 



where 



N B (w,P) 



1 



(A.4) 



(A.5) 



eP w - 1 e Py/p 2 +m 2 _ i ' 

As long asm^O the integrals ()A.3|) and (jA.4j) are not infrared divergent. Thus, before considering 
the zero mass limit an infrared cut-off p > must be introduced. In the limit m — * 0, Nb(w, (3) —> 
Nb(p; 13) and we obtain, 



7J (+) (x-?/;m,/5) m ^ = D { +\x-y) + — / — N B (p; (3) cosp{x + -y + ) + cosp{x~ -y~) , (A.6) 

27T p L J 

which is the two-point function of the free massless thermofield given by Eq. (j3.44jl . By consid- 
ering the free massless scalar thermofield theory as the zero mass limit of the free massive scalar 
thermofield theory, the infrared regulator p of the zero temperature two-point function should be 
identified with the infrared cut-off p' of the temperature dependent contribution in (|3.44jl . 
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